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ABSTRACT 



A squall line is investigated during a 3-hour period of its mature 
life. Upper air data from the mesonetwork in Oklahoma at 1700, 1830, 
and 2000 GMT is used to determine vorticity, divergence, and vertical 
motion at 50 mb intervals. These and other fields are computed and 
investigated in detail, and possible problems associated with them 
are discussed. Two methods of computing vertical motion are described 
in detail. 

A mesoscale wave with a wavelength of approximately 90 n.m. is 
found at 700 mb. It moves faster than the wind at 700 mb but slower 
than the squall line. This wave is investigated in detail and is 
hypothesized to be the cause and the terminator of the squall line. 

The kinematic vertical motion technique presented appears to 
give realistic fields. Appendix B describes the problems found in 
the vorticity method and how it can be used on the synoptic scale. 
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I. INTRODUCTION 



The severe convective storm, which produces death and destruction 
in the form of tornadoes, hail, and violent winds, has plagued man 
since the beginning of time. With the advance of the science of 
meteorology in the last century, man has attempted to increase his 
knowledge of these interesting yet deadly storms. House (1963) has 
stated that attempts were being made as early as the 1880's to anti- 
cipate the development of such storms. However, most of the work on 
the possible causes and structure of these phenomena has been done, as 
has most of the advances in meteorology, in the last twenty-five years. 

To date a considerable gap remains in man's understanding of the 
severe storm. A major problem has been his inability to measure the 
atmospheric parameters such as wind and temperature on the scale 
necessary to represent the severe storm. Within the last decade 
significant progress has been made with the introduction of obser- 
vations on the scale of these phenomena. This is the scale of motion 
(10 to 500 km) in which convective storms appear. 

The studies that have been made to date are at times contradictory, 
yet they all add to our general knowledge. Carlson and Ludlum (1968) 
point out that either a band of strong wind in the middle troposphere 
or strong vertical wind shear is necessary for the development of 
severe convective storms. Also, Koscielski (1965) found that the 
lack of an upper level jet was generally associated with no activity 
or less intense activity than that associated with such a jet. How- 
ever, Byers and Braham (1949) in the report of the Thunderstorm Project 
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in Florida and lately Asai (1964) have determined that the presence 
of vertical wind shear acts as a deterrent to convective growth of 
the ”air mass" variety. 

Newton (1966) has noted this different effect on severe storms 
and "air mass" cells, and proposes that the conflicting area is the 
intermediate zone where the two possibly intermingle or overlap. 

Takedo (1965) studied downdrafts and updrafts in storms. He 
noted that vertical wind shear seems to organize the release of 
energy in a convective system. When significant vertical wind 
shear was present, stronger upcurrents formed in front of the con- 
vective region than when they were absent. He indicates that vertical 
wind shear is very important for continuation of the system if the 
lower layer of the troposphere is very stable. Therefore, he con- 
cludes that vertical wind shear is necessary for continuance of a 
squall-type system. 

Newton (1966) and Alberty (1969) both note that the disposal of 
warm air in upper levels through the action of anvil blow-off, caused 
by vertical shear, is an important revitalizing mechanism of the 
severe storm. Alberty (1969) points out that although the strong 
vertical shear would be expected to force the convective system to 
lean in a down shear direction, the strong shear simultaneously yields 
an effect that enhances vertical motion near the strong updraft core. 
This increased vertical motion tends to overpower any leaning in the 
down shear direction since the environmental flow has less time to 
impart its horizontal momentum upon the vertically moving air parcels. 

Senshu (1961), in a study of a mesoscale disturbance in western 
Japan in 1955, proposed that the initial activity is caused by strong 
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convergence in the lower levels causing an upward current. Condensation 
aloft causes precipitation which falls, because of the rapidly moving 
system, sliglitly beiiind the convergence maximum. Some rain evaporates 
as it falls, cooling the air and causing the downdraft and a mesohigh 
to form. This increases the wind, and the local time change of wind 
velocity balances with the large pressure gradients preceding the 
system. This wind system then maintains the strong convergence in 
front of the high cell, thus maintaining the moving system. This 
mesoscale disturbance studied by Senshu (1961) moved at a speed of 
45 kts . 

Much of the work done in the United States on mesoscale systems 
can be attributed to T. Fujita of the University of Chicago. Fujita 
(1955) made a detailed study of two systems of squall lines in the 
central United States. Since only synoptic data was available in the 
area, Fujita developed a technique for converting time cross sections 
to space cross sections. By assuming that the local change of an 
element is small over an interval of an hour (using barograph traces, 
etc.), time sections were converted to space sections. This technique 
made available a much more detailed data network over the area covered 
by only synoptic-scale reports. Only surface data is used by Fujita. 

The assumption that the local change of an element is negligible over 
an hourly interval is reasonable in the middle and upper troposphere 
for synoptic-scale motions, but it seems somewhat debatable near the 
surface for sub-synoptic motions. Nevertheless, the scheme used by 
Fujita would seem to be a relatively good approach in areas where no 
other data is available. This directly points out the need for meso- 
scale networks, as noted by Fujita (1955, 1957). Fujita (1955) found 
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large surface divergences of 10 to 60 x 10 sec over areas from 
100 to 10,000 square miles in the vicinity of the squalls. He found 
that the storms develop as a small mesohigh and dissipate as a meso- 
depression after going through five stages. The movement of the 
squall lines was associated with a **pressure-surge** line, caused by 
the mesohigh. However, in one case, he noted that the system was 
initiated before a pressure-surge line was formed, suggesting the 
futility of seeking a **trigger*' mechanism in the surface pressure 
field at the initial stage. A similar study by Fujita and Brown 
(1958) again noted the mesosystem and its associated convective 
cells and stressed the above technique for tracking it (when a 
mesonetwork is not available), fortified by a dense raingauge net- 
work and continuous radar pictures. Fujita (1957) noted the direct 
relationship between mesosystems of 10 to 100 miles in size to severe 
storms, and proposed that similar features are found in other parts 
of the world, particularly in Japan. 

Considerable work has been done in Japan in recent years on 
mesoscale disturbances, some of which has already been mentioned. 
Matsumoto (1967) and Matsumoto et al (1967a, b) have worked with 
data obtained from the surface mesonetwork set up in Japan for the 
Japanese Heavy Snow Storm Project. The results of these studies will 
be noted in more detail later, but they indicated that squall lines 
and severe convective activity in general are initiated and controlled 
by the surface area of convergence. This same conclusion is found by 
Doroshenko (1960). However, he indicates that this convergence may be 
caused by a trough or other causes; therefore, it is still unclear if 
this surface convergence is the cause of the squall formation, only one 
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of a scries of events triggered by the cause of the squall line, or 
possibly t'ven a direct effect of tlie squall line itself. 

The possibility of quas i -s ta t ionary mesoscale waves of wavelength 
500 to 600 km setting off squall activity in the central United States 
has been proposed by Dirks, Mahlman, and Reiter (1967), They assumed 
constant static stability and vertical wind shear over the area and 
found that a stationary wave of varying wavelength (dependent on 
several factors) could be postulated which would account for the 
seasonal maximum of squall line formations which are initiated three- 
fourths of a wavelength downstream from the ridge over the Rocky 
Mountains, Considerable use was made of satellite pictures in this 
study to fill in the gaps left by the synoptic upper air network and 
special pilot balloon observations in the area. A major recommendation 
of this study was again the call for an operational upper air mesonet- 
work so that the possibility of such mesoscale waves might further be 
investigated. Because of the great expense of such a network, very 
few will likely be made available in the United States or elsewhere 
for years to come. 

An experiment by Tepper (1959) showed the importance of mesoscale 
data. Using nine experienced Weather Bureau forecasters, three were 
given good surface mesoanalyses and the other six provided the regular 
analyses available at regional forecast centers. The three given 
mesoanalyses completely outperformed the other six in all phases of 
the forecast. So, it was evident that a definite increase in fore- 
casting ability and understanding the weather would result if a 
sufficiently dense and financially practical surface mesonetwork 
could be employed. Such a network (called a Beta network) was set 
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up in central Oklahoma in the spring of 1960 by the National Severe 
Storms Project (NSSP)* It was used in the spring of 1960 and 1961 
by the NSSP staff (1961) to gather data on severe storms and their 
environment. This Beta network contained 36 stations 10 to 15 miles 
apart. A similar network was used in Japan in conjunction with the 
Japanese Heavy Snow Storm Project previously mentioned. Since no 
upper air mesonetwork existed in the area, aircraft were used in an 
effort to determine the middle and upper features of the severe weather 
which occurred. Computations from the data recorded by the aircraft 
indicated an area of convergence of magnitude 15 x 10 ^ sec ^ near 
800 mb in the vicinity of the activity. Again, because of the lack 
of conventional upper air data throughout the area and the problems 
associated with aircraft reconnaissance, it was strongly recommended 
that a mesonetwork of upper air stations be formed in the area. 

Because 1) central Oklahoma is a hotbed for convective activity, 

2) the Beta surface network was already established there, and 3) 
of the previously stated numerous appeals for an upper air mesonetwork 
somewhere in the United States, the National Severe Storms Laboratory 
(NSSL) of ESSA initiated a program of mesoscale upper air soundings 
in 1966 covering the south central Oklahoma area. This network con- 
sisted of nine upper air stations operating in the spring and early 
summer (maximum severe weather) which were to take soundings every 
one-and-a-ha If hours beginning at 1700 GMT on days when severe weather 
was expected (see Fig. 1). These stations are approximately 40 to 50 
n.m. apart. It was hoped that this data, combined with the Beta surface 
network and radar pictures from the WSR-57 radar located at Norman, 
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Oklahoma, would open new doors in the study of severe storms and their 
interactions with the environment. 

Considerable use of the data obtained on 28 May 1967 by this net- 
work (a case of isolated activity) has been made, most notably by 
J. C. Fankhauser (1968, 1969) of the National Center for Atmospheric 
Research (NCAR) . He investigated the thunderstorm which occurred 
that day using five consecutive sets of soundings covering the life- 
time of the storm, and described it using divergence, vertical motion, 
and other computed parameters. He finds in general that the resultant 
picture of the storm area using the above parameters shows strong sur- 
face convergence and upward motion aloft (average of 72 mb hr ^ in 
the vicinity of the storm). His figures indicate a strong convergence 
region near 500 mb which is of great interest, but he does not postu- 
late any ideas as to the cause or significance of this occurrence. 

His research shows that although the sounding system alone is not 
complete enough to make possible an intense study of singular thunder- 
storm dynamics, it does provide a reliable framework within which the 
possible cause and effect of thunderstorm systems can be studied. 

With the above conclusion in mind, this study was initiated using 
the data obtained by the NSSL mesonetwork on 30 May 1967. On this day 
an active squall line and numerous other cells occurred within the area 
In order to obtain as much detail as possible, a 13 x 11 grid was used 
for obtaining data used in computations. This grid covers the area 
from 96^45* to 99^45* West and from 33^45* to 36^15* North, cutting 
grid sections of 15* of longitude and 15* of latitude, and encom- 
passing the mesonetwork. Fig, 1 shows the majority of the grid area 
(the boundary points are dropped as they are lost in most of the 
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d i ffc'i'enc i ng scIumhos used in making compu ta t i ons ) . Wichita Falls (SPS) 
is noted in Fig. 1 but not noted in proceeding figures; however, SPS 
data is used at all times unless otherwise noted. Grid points (I, J) 
are referenced in later sections and are shown in part in Fig. 1. 

Tlie objectives of this study are as follows: 

(1) describe the convective system and its relation to the 
niesoscale environment, 

(2) develop a method for obtaining the vertical motion field, 

(3) seek an explanation for the squall line which occurs during 
the period. 

It is necessary to be able to fulfill the above requirements in 
t(iis and similar cases in order to understand the development and life 
cycle of destructive severe storms. 
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II. GENERAL SYNOPTIC SITUATION AND INDICATION 
OF SEVERE WEATHER PRODUCTION 



Since more information has been compiled on the synoptic scale con- 
cerning conditions necessary for severe weather production than on the 
mesoscale, it seems appropriate to look at the general synoptic situ- 
ation present in the central United States on 30 May 1967. The times 
of interest here are centered around 1830 GMT (1230 CST) , but since 
complete soundings on the synoptic scale are taken only twice daily 
at 0000 GMT and 1200 GMT, the general picture at 1200 GMT will pro- 
vide a preview of the conditions associated with the mesosystem which 
developed later in the mesonetwork area. 

Synoptic conditions necessary for the occurrence of severe weather 
have been noted by many authors. However, these vary and many are 
only general in nature. The Weather Bureau is presently in the process 
of publishing a new guideline for severe weather criteria. The general 
criteria used by the Weather Bureau were noted in an unpublished study 
by Koscielski (1965). He found that severe weather generally occurs 
when the following parameters are relatively strong: instability, 

500 mb positive vorticity advection, upper jet, low level moisture, 
and 850 mb temperature and moisture advection, Endlich and Mancuso 
(1968), using a computer objective analysis technique, have proposed 
that a combination of the above synoptic indicators and numerous 
kinematic quantities available through computer use by dynamical 
techniques will improve synoptic severe weather forecasting. But 
much work has yet to be done here. Probably the simpliest and best 
summary of severe weather indicators to date is that made by Dirks 
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(1969). He has noted that synoptic-scale upper level divergence and 
lower level convergence are sometimes associated with severe weatiier, 
however, such occurrences are not frequent enough to account for the 
general formation of mesoscale convective systems. He lists the 
following as conditions for the production of severe weather: 

(1) a warm, moist tongue at low levels; 

(2) general convective instability (S.I. < 4); 

(3) possible presence of an upper level short wave; 

(4) a low and/or upper level jet; 

(5) cold air advection aloft; 

(6) mid- tropospheric stable layer; 

(7) a cold or "dry front" is often present. 

Fig. 2 indicates the synoptic conditions prevalent in the central 
United States in the early morning hours on 30 May 1967. Fig. 2a 
shows that a cold front dipped southward into Oklahoma and Texas from 
a wave in east central Kansas and then trailed off northwestward as a 
stationary front. A "dry front", or Marfa front as it is often called, 
was located in far west Texas and eastern New Mexico. Southerly flow 
is indicated across east and central Texas converging into the meso- 
network area. A dying squall line is present just southeast of the 
mesonetwork over north central Texas, indicating the instability over 
the area of the previous night and evening. Fig. 2b shows the sur- 
face conditions at 1800 GMT, Fig. 2c shows the frontal configuration 
at 850 mb. Considerable thermal gradient is indicated behind the front 
in central Kansas. A warm, moist tongue is evident from west and south 
Texas across the mesonetwork to Kansas. The small cool pocket over 
north central Texas appears a direct result of the decaying squall line 
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FIG. 2a. Surface analysis for 1200 (KT on 30 May 1967 with 
fronts and instability lines. Mesonetwork is small square 
in Oklahoma. 
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FIG. 2b. Surface analysis for 1800 GMT on 30 May 1967 with 
fronts and instability lines. 



27 




FIG. 2c. 850 mb height and temperature analysis for 1200 G^f^ 

on 30 May 1967 with fronts and dewpoints (°C). 
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FIG. 2d. 500 mb height and temperature analysis for 1200 GMT 

on 30 May 1967 with jet axis. 
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FIG. 2e . 200 mb height analysis for 1200 GMT on 30 May 1967 

with jet axis. 
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thoro. Also, both Amarillo and Midland, reporting 30 and 25 knot winds 
ros pec t i VO ly , indicate a moderate wind maximum over the eastern Texas 
Panhandle and possible western Oklahoma, between the depression in 
eastern New Mexico and the warm ridge to the east. 

Fig. 2d shows a short wave trough at 500 mb with a wavelength of 
approximately 600 n.m. extending from northwest Kansas to northeast 
Texas. The double shafted arrow extending northeastward over the high 
plains of Texas represents the jet stream. Also notice the cold 
advection occurring over Oklahoma. 

Finally, Fig. 2e shows the same short wave trough at 200 mb 
located over Texas, indicating a well developed trough. The double 
shafted arrow again represents the jet stream with an apparent moderate 
maximum of 65 or 70 knots. 

Cold advection at 500 mb (Fig. 2d) over warm advection at 850 mb 
(Fig. 2c) indicates decreasing stability. To gain an understanding 
of whether or not this is significant, write the first law of thermo- 
dynamics in the form 



1 ^ ^ 

c dt dt c ’ 

P P 



Expanding and rearranging the above gives 



bT _ _1_ ^ ^ 9^ 
St c dt c 
P P 



\V*V T 
P 



Od 



Sp’ 



Taking of the above and interchanging the derivatives on the left 
side gives 



St VSzy 



Sz vc dty Sz vc y 
p p 



^ - 1 ? 1 ^' 



(1) (2) (3) 



(4) (5) 
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To determine the effects of advection on stability, analyze how 



(4) affects (1). For simplicity, measure \V*V^T along a streamline 
(natural coordinates) over the mesonetwork area at 850 mb (Fig. 2c) 
and subtract the value found at 500 mb (Fig. 2d) to get an estimate 
of the rate of change of stability over the mesonetwork area previous 
to the outbreak of squall line activity. 

Fig. 2c indicates a maximum of 10 knots diagonally across the area 
and a temperature gradient of -2°C/150 n.m. Fig. 2d indicates a wind 
of 30 knots and a temperature gradient of +1.5^C/l00 n.m. Using these 
values one obtains a rate of change of stability (or an increase in 



mb. OKC (only synoptic station in mesonetwork area on Figs. 2c and 2d) 
shows approximately a standard atmospheric lapse rate at 1200 GMT 
(0600 GST), Therefore, as the day progresses, slow destabilization 
takes place on the synoptic scale over the area, due to differential 
temperature advection alone. This would indicate a probable drop of 
-0.5^C km ^ hr ^ or slightly more as low level solar heating adds to 
the advective destabilization. However, this destabilization on the 
synoptic scale seems rather weak and, by itself, would seem to indicate 
possibly only weak showers. 

From Figs. 2a, 2c, 2d, and 2e it is apparent that most of the 
conditions necessary for severe weather development indicated pre- 
viously are present over the mesonetwork area at 1200 GMT. By 1600 
GMT, a squall line was developing over the western part of the meso- 
network, and considerable instability was present on the synoptic 
scale by 1800 GMT as is shown in Fig. 2b. 




over the layer 850 to 500 
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III. DATA AND ANALYSIS PROBLEMS 



Data lias been taken for three time periods (1700, 1830, and 2000 
GMT) from the NSSL mesonetwork in Oklahoma for 30 May 1967. The data 
was plotted (wind only) and analyzed (both for streamlines and isotachs) 
for 17 levels in the vertical (surface to 150 mb at 50 mb intervals). 
Only the winds are used for calculations in this study. Clark (1969) 
and Van Sickle (1969) have indicated that wind information is the most 
accurate data returned by the radiosonde for use in mesoscale studies. 
These winds have been smoothed by a (1, 2, 1) process over a 600 meter 
vertical layer in order to make them even more representative. Analysis 
on the mesoscale presents some problems which must be tackled with 
other than conventional synoptic thinking, as has been shown in detail 
by Fujita (1963), Basically, of course, time and space continuity is 
maintained to a good extent; however, at 50 mb intervals the exact 
analysis of the transition zones between low level systems and the 
general upper level jet flow are open to question and any of several 
answers may apply. However, even though analysis might be easier, 
a true picture of the atmospheric motion present during the period, 
complex as it is, would probably be obscured if more than a 50 mb 
increment was used. Time continuity was almost impossible in the low 
levels, even though only a seemingly small 90 minute interval occurs 
between soundings. It appears that mesoscale features in the low 
levels can change considerably within this time period, as indicated 
by Fujita (1963). 
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Surface data used for the three hour period (1700 - 2000 GMT) was 
obtained from the Beta network explained in Section I* Surface pres- 
sure tendency and rainfall were analyzed for 1-1/2 hour intervals. All 
radar presentations used in the study are from the WSR-57 radar at 
Norman, Oklahoma. All presentations utilized are at 0° elevation. 
Unfortunately, the Sensitivity Time Control (STC) component was in 
use at all times because of ground return, so no return is available 
within 25 n.m. of Norman. This is a drawback, especially at 1830 GMT, 
when it is apparent that cells are present in this area. 

All of the preceding analyses were done using the data from the 
NSSL computer output. The data was plotted on station and considered 
to be representative of the time of the sounding. As noted by 
Fankhauser (1969), the displacement of the balloon with altitude 
from the launch station, time lag of the different pressure levels 
from launch time, varying launch times, and differing ascent rates 
may all present problems on the mesoscale. A detailed investigation 
into these four possible problems showed that the latter two were not 
fundamental difficulties because of the close adherence to the rules 
set up for this mesoscale network. The former two were found to pro- 
duce problems only under a very few critical situations. Appendix A 
describes in detail the investigation made. In general, it was deter- 
mined that the data could be used in the conventional, synoptic manner, 
but corrections were made in the few critical situations mentioned. 
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IV. DltlSCRIPTION OF VORTTCITY, DIVKRGKNCK, AND 
VERTICAL MOTION CALCUI^TIONS 



A. VORTICITY AND DIVERGENCE 

In ordc'r Lo host rc‘prc‘sc'nt t!ie mososcalc' (‘nvironment in soutf)- 
WGStorn Oklahoma on 30 May 1967, the patterns of relativ(‘ vorticity 
and divergence, which are well understood and described on the synoptic 
scale but much less understood or described on the mesoscale, were 
deemed necessary quantities. 

The definition of relative vorticity is 



" = ^ ^ 

^ r hx dy^ 

with all indicated differentiations to be performed on a surface of 
constant pressure. 

Relative vorticity was determined using a centered difference 
scheme. Non-centered differencing was used at the boundaries. 
Absolute vorticity was also calculated. 

The definition of divergence is 



Both relative 
the vertical by a 
Laplacian .15-. 15- 



DIV - V-'.V • + #;. 

OX oy 

vorticity and divergence fields were smoothed 
.2-. 6-. 2 smoother, and in the horizontal by a 
.4-.15-.15 smoother. 



i n 



B. DESCRIPTIONS OF POSSIBLE VERTICAL MOTION SCHEMES 

The determination of reliable vertical motions in the atmosphere 
is of fundamental importance in a number of meteorological problems. 
Because there is no generally accurate method of computing vertical 
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motion, no method can be tested completely. However, qualitative 
methods such as comparing computed vertical motions to radar pictures, 
rainfall records, and local pressure change fields can give acceptable 
verification to a system. Also, such principles as conservation of 
potential vorticity can be used for verification purposes. It there- 
fore is possible to verify areas of upward motion where precipitation 
occurs to some degree of success , but areas of downward motion are 
difficult if not impossible to substantiate. However, the upward 
motion areas are the weather producing areas, and therefore one can 
conscientiously concentrate on obtaining accuracy in these regions. 

Many techniques and methods have been developed and used on the 
synoptic scale. Ballenzweig (1955) reviewed many of these, noting 
that they can generally be broken down into three categories: 
adiabatic, kinematic, and vorticity methods. All have their dis- 
advantages and advantages, but particular ones seem more suited for 
use on the mesoscale. 

Panofsky et al (1959) compared several methods of determining 
vertical motions on the synoptic scale and concluded that the kine- 
matic and omega equation techniques were best. However, a necessary 
condition for the validity of any omega equation method is that the 
balance equation be valid, and it is only valid if the Rossby number 

(R^) is less than 1. Using appropriate mesoscale motions such as 

-1 -4 -1 5 

V 20 m sec , f ~ 10 sec , and L (wavelength) ^ 100 km = 10 

meters, one sees that 




Therefore, the balance equation is not valid, so no omega equation 
method is valid on the mesoscale. 
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Tho use of an adiabatic or thermal equation technique would appear 
less desirable because 1) of the frequently observed weak thermal 
gradients on tlie mesoscale, 2) most of these techniques give time 
averaged vertical motions whereas instantaneous values would be of 
greater interest for mesoscale studies, and 3) in general, the mass 
field is difficult to measure accurately on this scale of motion. 

Van Sickle (1969) used this type of technique and noted that the 
problems associated with such an approach made it undesirable. 

The above critiques lead to an init ial choice of a kinematic 
type technique. But it is desirable to have two methods for com- 
puting vertical motion for cross checking. The previously indicated 
reliability of wind data and the inapplicability of other techniques 
lead to a choice of a vorticity method. Although the twisting term 
is generally considered negligible on the synoptic scale, it apparently 
must be retained on the mesoscale (Matsumoto et al (1967a)), so the 
entire vorticity equation was chosen. A scheme similar to the one 
used here was used by Yanai and Nitta (1967) in a synoptic study of 
a Caribbean easterly wave. They retained the troublesome twisting 
term, as it was found not to be negligible in the upper troposphere. 
They obtained a convergent scheme, which gave hope that the important 
twisting term could also be included successfully on the mesoscale. 

C. THE ADJUSTED KINEMATIC TECHNIQUE 

The major problem with the kinematic method is that errors are 
additive and as successive omegas are determined, as one integrates 
from the surface up, the method tends to "blow up" near the top of 
the atmosphere. The scheme used here, similar to the one suggested by 
Mahlman (1968), is to force the average divergence between two levels 



37 



(generally a 50 mb increment except for the first layer - surface to 
900 mb) to conform to the vertically averaged divergence determined 



previously by an integral constraint. This must be the case if the 
computed vertical motions are to agree with the upper and lower 
boundary conditions. 

The kinematic technique used here was developed in the following 
manner. First, the divergence is found at all interior points of the 
grid area and smoothed. Then the average divergence for each vertical 
column is found by mass -weight ing the 17 levels of divergence and 
averaging . 

Integration of the continuity equation between the surface and 
150 mb gives 



where sfc = surface and D = vertically averaged isobaric divergence. 

Surface pressures were determined from topography and the standard 
atmosphere, and surface winds were taken to be those at 950 mb. This 
was done in an attempt to get a representative surface omega at any 
time . 

Assuming that the vertical motion at 150 mb is zero, or approxi- 
mately so (this seems a good assumption since examination of all 
soundings indicated that the 150 mb level was above or at the tropo- 




( 1 ) 



vw 



w 

pause for all cases), one can rewrite (1) solving for D as 




( 2 ) 



where is determined by the terrain effect, or 




38 



This surface omega is used in both tlie adjusted kinematic and complete 
VO r I icily methods. 

Equation (2) then defines the vertically averaged divergence for a 
column of air to be that determined by the surface omega and the depth 
of tlie column in millibars. One may average all the divergences com- 
puted over the column, but this average must equal that value determined 
bv (2). Therefore, (2) acts as a ’’constraint" on the divergence fields, 
and will hereafter be referred to as the integral constraint. It is 
well known that divergence in itself is difficult to obtain accurately 
on the synoptic scale, since and ^ are often nearly equal and of 
opposite sign. Therefore, a slight error in either term can cause a 
relatively large error in the divergence. This problem may not be as 
great on the mesoscale, since average divergences appear to be an order 
of magnitude larger than those on the synoptic scale (noted earlier), 
but it still exists. Therefore, since the integral constraint is 

determined from accurately measured quantities, forcing individual 

w 

divergences (through the use of D) to meet this criteria by correcting 
them in some manner would seem to give divergence fields with less 
error that more accurately represent the state of the atmosphere. 

Reiter (1963) and Clark (1969) have shown that radiosonde errors 
tend to increase linearly with height. This assertion indicates that 
vertical velocities determined from the computed divergence patterns 
decrease in accuracy as one proceeds upward, so a normalized scheme 
for a smaller correction to the lower level divergences and a greater 
correction to the upper fields was developed. This method is similar 
to the one used by Fankhauser (1969) which was developed by O’Brien 
(1969). Since 17 levels have been used here, a mid-level (level 9 
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or 550 mb) was chosen as the level to which the total correction was 
applied to the divergence field. This total correction (CORK) is the 
difference between the integral constraint and the vertical average of 
the mass-weighted divergences of all 17 levels. By using a vertical 
indicator K, the divergence fields in the vertical were corrected by 
the following: 

DIV (NEW) = DIV (OLD) + | (CORK) 

where K ranges from 1 at the surface to 17 at 150 mb, 

and CORK = D - computed vertically averaged divergence. 

(NOTE: Abbreviations such as CORK are the ones used in the computer 

program at the end of this text and are used here to maintain 
continuity, ) 

The final corrected divergence fields were then used to calculate 
the omega fields at all levels. 

Calculation of vertical motions by the adjusted kinematic technique 
was done in five steps in order to compare the preference of one answer 
to the other. These steps were: 

(1) no corrections to divergence (true kinematic); 

(2) total correction; 

(3) K/ 9 times total correction; 

(4) K/9 times total correction with horizontal divergence smoother; 

(5) K/9 times total correction with horizontal and vertical 
divergence smoothers. 

The above steps correspond to Figs. 3-7 for time 1700 GMT at 
point (5, 5) (the Wichita Mountains). This was the vicinity of the 
squall line at 1700 GMT. Examination of Figs. 3-7 shows that a 
general improvement in vertical consistency occurs as one progresses 
through the figures. The extremely large values of divergence near 
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FIG, 3, Vertical profiles of divergence ( ) and vertical 

motion ( ), determined by the true kinematic method with 

no corrections, for grid point (5, 5) (Wichita Mountains) at 
1700 GMT on 30 May 1967. 
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FIG. 4. Vertical profiles of divergence ( ) and vertical 

motion ( ), determined by the completely corrected kine- 
matic method, for grid point (5, 5) (Wichita Mountains) at 
1700 GMT on 30 May 1967. 
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FIG. 5. Vertical profiles of divergence ( ) and vertical 

motion ( ), determined by K/9 times the complete correction 

for the kinematic method, for grid point (5, 5) (Wichita 
Mountains) at 1700 GMT on 30 May 1967. 
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FIG. 6. Vertical profiles of divergence ( ) and vertical 

motion ( ), determined by K/9 times the complete correction 

for ‘ the kinematic method with a horizontal divergence smoother, 
for grid point (5, 5) (Wichita Mountains), at 1700 GMT on 
30 May 1967. 
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FIG. 7. Vertical profiles of divergence ( ) and vertical 

motion ( ), determined by K/9 times the complete correction 

for the kinematic method with a horizontal and vertical diver- 
gence smoother, for grid point (5, 5) (Wichita Mountains) at 
1700 GMT on 30 May 1967. This is the adjusted kinematic 
technique. 
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200 mb are adjusted to more realistic values, but the general patterns 
of divergence and vertical motion are unchanged. Therefore, it was 
concluded that vertical motions determined by step (5) above most 
closely represent the real atmosphere. All vertical motion fields 
presented in this study were determined in this manner. 

It is interesting to note the area of maximum convergence which 
appears near the 600 mb level in Figs. 3-7, which appears to be a 
prevalent condition in all the shower areas as noted in a later section. 



D. THE COMPLETE VORTICITY METHOD 

The vorticity equation in pressure coordinates (neglecting friction) 
is 

dC, 



a . V 7 u/ a. ^ ^ ^ _ 



The following derivation gives the basis for the vorticity 
method chosen. Expansion of (3) gives 






Sc, 



SbD Sv cb 3u 



V. • + + CD V Q V-W + 3 — - 3 — 3 ;— = 0 

ot ^a op a ox op oy op 



SCa Sc^ 

Since , (4) becomes 

^^r , , ^^a , r r-r u/ 1 ^ ^ cb 5 u _ 

■^7 ^ W-VC + CD + = 0 

ot ^a Op a ox Op c3y c3p 

CtD 

By making the continuity equation substitution of V*\V = - 
(5) becomes 



Sc 



Sc 



Sp Sp 

Now, combining the bracketed terms in ( 6 ) gives 

sc 



''^1^ , \\i jrrr , 2 S /O) 'N , ctJO Sv So) Su _ 



(3) 



(4) 



(5) 



... r ^a . Stu”! , Stu dv cto Su ^ 

^ + w-vc^ + ^ - C, 5 ;J + ^ ^ - - 5 ; ^ - 0 . ( 6 ) 



(7) 
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Bn' rt'arranging , one obtains 






^ CD 



St 



■ r ^ Sv Sd Su 



Sx Sp Sy Sp 

,2 ,2 



( 8 ) 



Integrating (8) gives 






O) 

u 

a 



r -z—^ + w-vr 




i'” ^ ^ ^ ^ 

5p + I 6p. 



(9) 



^a 

P P 



Since omega occurs on both sides of the equation, a first approxi- 
mation must be made and then convergence must occur for the equation 
to be satisfied. 

The method of approach used here is to define 

3"^ + W-VC 
FUNl = 2 



FUN2 



ch) ^ ^ ^ 

= “Sx 5p by Sp 

c,^ ■ 



For the first approximation, the same type scheme is used as was 
used in the adjusted kinematic method. Using only FUNl, find the 
average FUNl for a vertical column. Then determine the average FUNl 
which satisfies the integral constraint. This is done by substituting 
FUNl and FUN2 into (9) to get 

03 ^ 03t 



sfc 



c. 






sfc 



] icf) VAWW WWW 

= FUNl Ap + FUN2 Ap, 

^150 



( 10 ) 



47 



WAWV 

where FUNl = vertically averaged FUNl, 

WWW 

FUN2 = vertically averaged FUN2 , and 

= P ^ “ 150 mb . 

s fc 

As in the adjusted kinematic method, if 150 mb) is assumed to 

be zero, (10) can be rewritten as. 



CO 



c -pc WWW 

= FUNl + FUN2. 



s fc 



( 11 ) 



vww 

The left side of (11) then acts as a constraint on both FUNl and 

wwv 

FUN2 , and will hereafter be referred to as the integral constraint. 

Since (11) contains two unknowns, the following procedure was used in 

vwwv 

solving it. First, one may determine FUNl by the integral constraint 
WWW 

by assuming FUN2 is zero. Then, vertically average the FUNl values 

at all 17 levels through the column by mass -weight ing them. The 

wwv 

difference between this averaged FUNl and that determined by the 

integral constraint is called CORVFl (see Appendix B). This CORVFl 

then is the average correction which must be made to the column so 

that the integral constraint is satisfied. The same K/9 times CORVFl 

correction technique as was used in the adjusted kinematic technique 

is used to correct the FUNl values at each level. 

Using these corrected FUNl values, and assuming that all FUN2 

terms (the twisting term) are zero (as is done on the synoptic scale), 

one may use (9) to integrate upward from the surface a level at a time 

to get a first approximation for the omega field. 

The next step is to determine the fields of FUN2 , This is done 

by using the first approximation omega fields just found. In satis- 

www 

fying the integral constraint (11), it was determined that FUN2 must 
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be zero. The FUN2 values at all 17 levels through the column are 
vertically averaged using mass-weighting, as was done in FUNl. So, 

WWW 

to meet the integral constraint, this vertically averaged FUN2 must 
be zero. If it isn't, the average correction, called C0RVF2 , is used 
to correct the FUN2 fields in a manner similar to the one used to 
correct the FUNl fields. The difference is that neither the top nor 
bottom FUN2 fields are corrected. The lowest FUN2 field is specifically 
determined by the terrain effect, and the top FUN2 field (at 150 mb) 
was determined so that the omega field here was zero. Both of these 
then must be kept intact (see the Fortran IV computer program - sub- 
routine VOROMG for details). 

With the corrected FUNl and FUN2 fields (9) may be used to 

determine the omega or vertical motion field. This field must then 

be compared to the first approximation field obtained from the FUNl 

field alone. If any set of omegas in the two omega fields differ by 

some predetermined epsilon, then the last omega field replaces the 

wwv 

first approximation field and one recycles through the FUN2 correction 
procedure and determines a new omega field. This field is checked 
against the old one, and the procedure continues until a satisfactory 
omega field is found. 

The results of this vorticity method are found in Appendix B. 

The Fortran IV computer program used to determine all the vari- 
ables described in this section is found at the end of this text. 

A determined effort has been made to provide sufficient comment 
cards to make this program self explanatory. The adjusted kinematic 
technique is the subroutine KINOMG, and the vorticity method is the sub- 
routine VOROMG. The nomenclature used in the above explanations is 
that which is used in the computer program. 
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V. RESULTS 



A. POSSIBLE PROBLEMS 

Since the storm began near 1530 GMT and was still active as a 
squall line at 2000 GMT, even though beginning to dissipate, it is 
considered to be a mature, active squall line throughout the 1700 to 
2000 GMT period of study. The severity and intensity of the squall is 
indicated by the report of hail 20 miles southeast of Altus (LTS) at 
1600 GMT. Therefore, a close look at the mesoscale environment in and 
around the squall line area becomes a necessity. However, because the 
squall line has reached maturity before the mesonetwork is activated, 
it has definitely begun to interact with and possibly reshape the 
mesoscale environment in which it is embedded. The effect the squall 
line, or any such mesoscale phenomenon, has on its environment is yet 
a mystery. Data which was possibly under the effect of the squall 
line was used with this in mind, if used at all. The above should 
always be kept in mind when evaluating the resultant fields which 
describe the mesoscale environment. 

B. THE MESOSCALE ENVIRONMENT AT 1700 GMT 

At 1700 GMT, the active squall line was oriented generally north- 
south along a line just west of Wichita Falls (SPS) and Fort Sill 
(FSI) . Fig. 8 shows the streamline and isotach fields for repre- 
sentative levels at this time. General cyclonic curvature is found 
over the squall line in the low levels and anti-cyclonic curvature 
in the middle and upper troposphere. A wind maximum at 200 mb over- 
laps a minimum at 500 mb over the squall line indicating large shears. 
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FIG. 8a. Streamlines ( — >) , isotachs ( ) (m sec" ), 

and the axis of the squall line (***) for 700 mb at 1700 
GMT on 30 May 1967. 
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FIG. 8b. Streamlines ( — >) , isotachs ( ) (m sec ^) , 

and the axis of the squall line (***) for 500 mb at 
1700 GMT on 30 May 1967. 
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FIG. 8c. Streamlines ( — ^>) , isotachs ( ) (m sec ), 

and the axis of the squall line (***) for 200 mb at 
1700 GMT on 30 May 1967. 
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Fig. 9 shows the squall line and the lower tropospheric fields associ- 
ated with it. Convergence at 950 mb is indicated over most of the 
mesonetwork area, with the maximum area of convergence oriented parallel 
to and just east or downwind of the squall line. Although this area of 
convergence is similar to that noted by Matsumoto (1967) and Matsumoto 

et al (1967a, b), it is not so strong (only a small area of magnitude 

-4 -1 

10 sec ) and generally precedes the squall line rather than encom- 
passing it. This seems plausible since the squall line is moving 
rather rapidly eastward, as is noted later in this section. The 
vertical motions shown here, as those shown for all times, are those 
determined by the adjusted kinematic technique (as noted in Fig. 7). 

The vertical motion field at 800 mb seems rather complicated and 
difficult to explain. The relative humidity pattern shows that the 
most moist air (90^ or greater) at 850 mb is also oriented parallel 
to and just east of the squall line. A large gradient of humidity is 
noted westward of this area, with rather dry air (30^ or less) indi- 
cated behind the squall line near the western edge of the grid area. 
Therefore a good energy source in the form of low level, convergent 
moist air is apparently present to sustain the squall line. 

Fig. 10 shows the squall line and the mid-troposphere divergence 
and omega fields. Also noted is the 1-1/2 hour rainfall record. Of 
most importance here is the fact that all the convective activity is 
in the area of upward vertical motion at 450 mb. This level was 
chosen so that comparisons could be made with Fankhauser (1969). He 
found that all the convective activity studied in the area two days 
before was also determined by the vertical motion at 450 mb. He 
determined an average magnitude of near 75 mb hr ^ in the vicinity 
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FIG. 9. Squall line and low level fields at 1700 GMT 
on 30 May 1967. 

950 mb divergence ( ) 10 ^ sec ^ ^ 

800 mb vertical motion ( ) mb hr 

850 mb relative humidity (-•-) ^ 

700 mb trough (==) 

Radar return (■«— ) 
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FIG. 10. Squall line and mid-troposphere fields at 
1700 GMT on 30 May 1967. 

550 mb divergence ( ) 10 ^ sec 

450 mb vertical motion ( ) mb hr 

3-hour rainfall pattern centered at 
1700 GMT (-•-) hundredths of inch 
700 mb trough (==) 

Radar return (a™*) 



56 



of tlie showers. Vertical motions in Fig. 10 appear somewhat greater. 

The maximum of near 240 mb hr ^ west of SPS seems somewhat large, but 
reanalysis in this area produced no change. The vertical motion maxi- 
mum axis appears parallel to but slightly behind or upwind of the squall 
line. This is somewhat surprising, but a possible explanation will be 
made later. Areas of maximum rainfall are in good agreement with this 
vertical motion pattern. Strong convergence is found at 550 mb over 
and to the rear (upwind) of the squall line. Similar findings were 
noted by Fankhauser (1969) and the NSSP staff (1961), however no 
attempt has been made in either of these cases to determine the sig- 
nificance, if any, of this strong convergence pattern at this level. 

C. THE MESOSCALE ENVIRONMENT AT 1830 GMT 

By 1830 GMT, the squall line has remained intact and moved rapidly 
eastward to the vicinity of Tinker (TIK) and Chickasha (CHK) . Fig. 11 
shows the streamline and isotach pattern over the area at this time. 

At 700 mb and to some extent at 500 mb cyclonic curvature is found 
just west of the squall line and anticyclonic curvature over the squall 
line. At 200 mb, as at 1700 GMT, anticyclonic curvature is noted in 
the vicinity of the squall. Also, as at 1700 GMT, an area of maximum 
winds at 200 mb appears in the vicinity of a wind minimum at 500 mb 
over the squall line. 

Fig. 12 shows the squall line and low level fields at this time. 
Unfortunately, since the line is over the WSR-57 radar at Norman, the 
activity within 25 n.m. of Norman is deleted from the return by the 
STC component in operation. However, the rainfall record in Fig. 13 
indicated that there was activity within this area during the period. 
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FIG. 11a. Streamlines ( — >) , isotachs ( ) (m sec ), 

and the axis of the squall line (***) for 700 mb at 
1830 GMT on 30 May 1967. 
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FIG. 11b. Streamlines ( — ^>) , isotachs ( ) (m sec , 

and the axis of the squall line (***) for 500 mb at 
1830 GMT on 30 May 1967. 
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FIG. 11c. Streamlines ( — ^>) , isotachs ( ) (m sec ^), 

and the axis of the squall line (***) for 200 mb at 
1830 GMT on 30 May 1967. 
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FIG* 12. Squall line and low level fields at 1830 GMT 
on 30 May 1967. Note: No radar return OKC-TIK area 

because STC component was in operation. 

950 mb divergence ( ) 10 ^ sec ^ ^ 

800 mb vertical motion ( ) mb hr 

850 mb relative humidity (-•-) ^ 

700 mb trough (==) 

Radar return 
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FIG. 13. Squall line and mid-troposphere fields at 
1830 GMT on 30 May 1967. Note: No radar return 

OKC-TIK area because STC component was in operation. 

550 mb divergence ( ) 10 ^ sec 

450 mb vertical motion ( ) mb hr 

3-hour rainfall pattern centered at 
1830 GMT (-*“) hundredths of inch 
700 mb trough (==) 

Radar return («■!■*) 
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Fig. 12 shows a weak area of convergence in the vicinity of and just 
preceding the squall line. A strong area of divergence is located 
north of Fort Sill (FSI). This is directly associated with the strong 
downward vertical motion in the area, and appears to be well confirmed 
by the relative humidity, which is a minimum in this area. Again, high 
relative humidity appears to be associated with the area of convergence, 
but this is not as clear as was the case at 1700 GMT. The upward 
vertical motion maximum at 800 mb is in the area of the squall line. 

Fig. 13 shows the mid -tropospheric environment and squall line. As 
at 1700 GMT, all the echoes are within the upward motion area at 450 mb. 
The field seems very realistic here, with two maximum upward motion 
areas of 100 mb hr These don’t exactly overlay the two maximum 

rainfall areas, but they are connected by the same axis which connects 
the rainfall maximum. Since the vertical motions are instantaneous 
and the rainfall occurs over a 1-1/2-hour interval, this relationship 
seems very encouraging. 

As at 1700 ®1T, large convergences appear behind or upwind of the 
squall line at 550 mb. 

D. THE MESOSCALE ENVIRONMENT AT 2000 GMT 

By 2000 GMT the squall line is just passing out of the eastern 
end of the grid area. It is aligned generally north-south from 
96^30’W in the north to 97^00*W in the south. Fig. 14 shows the stream- 
line and isotach fields for representative levels at this time. The 
curvature at all three levels is slightly anticyc Ionic over the squall 
line at this time. A wind maximum remains at 200 mb over the far 
northern end of the squall line. At 500 mb a weak tongue of maxi- 
mum wind is indicated over the northern end of the squall line. At 
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FIG 14a. Streamlines ( — >) , isotachs ( ) (m sec 

and the axis of the squall line (*^) for 700 mb at 
2000 GMT on 30 May 1967. 
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FIG. 14b. Streamlines ( — ^>) , isotachs ( ) (m sec ^), 

and the axis of the squall line (***) for 500 mb at 
2000 GMT on 30 May 1967. 
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FIG. 14c. Streamlines ( — >) , isotachs ( ) (m sec ^), 

and the axis of the squall line (***) for 200 mb at 
2000 GMT on 30 May 1967. 
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700 mb strong cyclonic curvature or troughing is located to the west 
of the squall between CHK and TIK , with a wind maximum of 19 m sec ^ 
noted at CHK. 

Fig. 15 shows the radar return and the lower tropospheric fields 
at 2000 GMT. The squall line at the far eastern edge of the area is 
notably more diffuse and the individual cells it contains are smaller 
and weaker than at either of the previous times. However, note the 
"explosion" of minor, microscale echoes throughout most of the area 
at this time. No clearcut reason is evident for this random breakout, 

but it may be due to local heating effects since this is the maximum 

heating time of the day (1400 CST). Since calculations of vertical 
motion and divergence end at 97^00*W, only the southern end of the 
squall line is covered by these parameters. Fig, 15 shows strong 
sinking motions and divergence over most of the area in the lower 
levels. Weak upward motions and weak convergence is indicated in 
the far northwestern and southwestern parts of the area. Again, 
relative humidities decrease from values greater than 80^ in the 

far east to a large circular area of 30^ or less near the center of 

the area (FSI and CHK). 

Fig. 16 shows the mid-tropos pheric environment at 2000 GMT. Most 
of the area is covered by upward motion at 450 mb. However, an L- 
shaped area of sinking motion extends generally from the northern 
edge of the area south along 98°W to CHK then eastward. A maximum 
of 20 mb hr ^ sinking is noted just north of CHK. A few cells exist 
in the northern edge of the area where sinking is very weak, but as a 
rule the area of sinking motion is practically clear of echo returns, 
expecially when contrasted to adjacent areas to the south and northeast. 
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FIG. 15. Squall line and low level fields at 2000 GMT 
on 30 May 1967. 
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FIG. 16. Squall line and mid-troposphere fields at 
2000 GMT on 30 May 1967. 

550 mb divergence ( ) 10 ^ sec ^ ^ 

450 mb vertical motion ( ) mb hr 
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As has been noted at all three times and also by Fankhauser (1969), 

in parctically all cases echoes exist only if there is upward motion 

present in the middle troposphere. However, the echoes present may 

be very weak and shallow, regardless of the vertical motion at lower 

levels. Fig. 16 shows convergence over the whole area at 2000 GMT, 

-4 -1 

with values of order 10 sec in the area just behind the squall 
line. The rainfall for the period 1915 to 2045 GMT noted in Fig. 16 
indicates the majority of rainfall during the period was caused by 
the squall line (maximum of 0.50 in at the eastern edge). Therefore, 
it is impossible to correlate it with the vertical motion picture at 
2000 GMT, except possibly on the southern end of the squall line which 
is still in the area. Here a maximum area of upward motion of 40 mb hr ^ 
over Pauls Valley (PVY) is in good agreement with a bulge of the 0.25 in 
isopleth of rainfall. 

E. EVIDENCE OF A 700 MB MESOSCALE WAVE 

1 . As a Possible Cause of the Squall Line 

The large areas of convergence found in the middle troposphere 
and the apparent association of celular activity and the 450 mb upward 
motion area suggest that a possible mechanism for initiation and main- 
tenance of the squall line system might be found in the middle troposphere. 

As mentioned previous ly, Mats umoto et al (1967a) found strong low 
level convergence definitely associated with organized convective 
activity, but they also noted pressure anomoly lines associated with 
such activity. They noted that these pressure anomoly lines moved in 
the direction of the 700 mb wind with a speed of about 45 knots. 

Matsumoto did not make streamline analyses above the surface, and no 
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attcMiipt to determine the possible presence of a 700 mb mesoscale wave 
was made. 

Investigation of the middle and upper troposphere streamline fields 
(every 50 mb) showed that a trough existed in the 700 to 600 mb layer 
at all three times to the rear of the squall line, and is best noted 
at 700 mb. It can be seen in Figs. 8a, 11a, and 14a. A definite 
trough was not found at any other levels aloft. The trough moved in 
the same direction (ENE) as the squall line but slowly lost ground to 
it. The wavelength of the trough appeared to be approximately 90 n.m. 

The trough moved at an apparent speed of 20 knots from 1700 GMT to 1830 
GMT, and 30 knots from 1830 GMT to 2000 GMT. Since the trough was 
moving generally from west to east, a comparison of its speed to the 
u component of the wind was necessary to determine its movement rela- 
tive to the air parcels at that level. The average u value for 700 mb 
over the entire grid for 1700 GMT was 11 knots, at 1830 GMT was 13 knots, 
and at 2000 GMT was 10 knots. Therefore, the trough moved at an aver- 
age speed of 25 knots from 1700 to 2000 GMT, but the air parcels at 
700 mb moved only at an average of 12 knots. The trough then moved 
at a phase speed over twice as fast as the air parcels. At the same 
time the squall line moved at 30 knots from 1700 to 1830 GMT and 37 knots 
from 1830 to 2000 GMT, or an average of 33 knots for the 3-hour period. 
Therefore, the squall line slowly moved away from the trough during this 
period. The rapid rate of speed of this squall is hard to correlate with 
any particular level, but seems to best fit the 450 to 500 mb level where 
the average eastward component of the wind for the 3-hour period was 
near 25 knots, a maximum for all layers. Since this is only the east- 
ward component of the wind (500 mb wind direction is generally from 
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the southwest) and the squall moved slightly ENE , it seems that this 
level is best correlated with the moving, mature squall line during 
this period. 

The trailing 700 mb trough moving slower than the squall but 
faster than the wind at 700 mb poses the interesting question of 
whether it is a possible cause or effect of the squall, A first 
indication that it is a cause rather than effect is the fact that its 
movement appears nearly independent of the movement of the squall, 
Another indication is the observation made earlier that the shower 
activity appears to correlate well with the vertical motion pattern 
in the middle troposphere. How does this 700 mb trough appear in 
other fields? Is it supported or not? That is to say, is it really 
there? 

It has been seen that the 700 mb trough appears in the streamline 
field. An analysis of the surface pressure tendency was done for 
1-1/2-hour periods centered on 1700, 1830, and 2000 GMT. These 
analyses are seen in Figs. 17a, 17b, and 17c respectively. These 
analyses showed a definite negative surface pressure tendency line. 
This line paralleled the 700 mb trough axis determined by the stream- 
lines and was positioned about 10 n,m, in front of it at 1700 and 
1830 GMT, At 2000 GMT the pressure tendency line was somewhat more 
complex and appeared to be misshapen by the random local cells in the 
area, but as a whole still preceded the 700 mb trough position deter- 
mined by the streamline field by an average of 10 to 15 n,m. This 
seems to support the proposed trough at 700 mb, and indicates that at 
least to this extent it is affecting conditions at the surface. 
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FIG. 17a, The 700 mb trough, squall line axis, and 
surface pressure tendency at 1700 GMT. 
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FIG. 17b. The 700 mb trough, squall line axis, and 
surface pressure tendency at 1830 GMT. 
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FIG. 17c. The 700 mb trough, squall line axis, and 
surface pressure tendency at 2000 GMT. 
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It has been shown that the trough is moving faster than the air 
parcels at 700 mb. Therefore, with respect to the trough, air parcels 
are moving toward and through it from the east, and away from it to 
the west, even though both the trough and air parcels are moving in 
the same absolute direction. Since the trough is generally not in a 
shower or echo area, it is reasonable to assume that adiabatic motion 
takes place, or that air parcels moving into the trough are travelling 
on an isentropic (0) surface. It is well known that potential vorticity 
(P) is conserved on a 0 surface, or 



dt 




= 0 . 



Expansion of the above expression gives 
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where C = V K - — , in natural coordinates. 

'' r s s on 

Inspection of Figs. 8a, 11a, and 14a shows that in the last term, as 
a parcel moves into the trough, the curvature part of the term is much 



greater than the shear part ^ ]• a parcel moves into 

the trough the curvature term increases (V^K^), or relative vorticity 
(C^) increases. As it moves westward out of the trough, the curva- 
ture term, which is a maximum in the trough, decreases, or decreases. 
Consequently, should be a maximum in the trough. 

For potential vorticity to be conserved, the above expansion shows 
that an increase in relative vorticity must be accompanied by a de- 
crease in stability An analysis of ( 7 ^- 7 ), which is a form 

of static stability proportional to (^- made for all three 
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times. If a relative vorticity maximum does occur along the proposed 
trough, then it should be accompanied by a stability minimum. 

Therefore, if the trough exists it should contain a minimum of 
stability in association with the maximum of relative vorticity along 
its axis. Figs. 18, 19, and 20 show the relative vorticity field, the 
static stability field, and the trough axis at 700 mb at 1700, 1830, 
and 2000 GMT respectively. At 1700 GMT the trough is near the edge 
of the grid area and is not too clearly defined, but at 1830 and 2000 
GMT it is unmistakably there. Therefore, the static stability mini- 
mum and vorticity maximum are in the trough to the rear of the squall 
line at all three times. 

Following the logic above, it is seen that if relative vorticity 
increases in front of the trough, the vorticity equation on a constant 
pressure surface (neglecting friction and assuming the twisting term 
is negligible) 



^ C + C V-W = 0 
dt ^a ^a 

shows that there must be an increase in convergence. This convergence 
would be most pronounced near the inflection point in front of the 
trough. This convergence would normally cause a downward motion at 
the inflection point and below it, and an upward motion above 700 mb. 
The reverse would be true behind the trough. At all three times the 
divergence fields at 700 mb show a weak area of convergence in front 
of the trough in the area of the inflection point. From 750 mb to 
the surface there is generally divergence, but this is not so clear 
near the lower levels because of the terrain effect. However, the 
area of convergence is stronger above 700 mb with a maximum near 550 
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FIG. 18. 700 mb trough, vorticity, and static stability 

at 1700 GMT on 30 May 1967. 
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FIG, 19. 700 mb trough, vorticity, and static stability 

at 1830 GMT on 30 May 1967. 
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FIG. 20. 700 mb trough, vorticity, and static stability 

at 2000 (KT on 30 May 1967. 
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mb. This is more true at 1830 GMT than at the other two times. This 



would naturally cause a larger depth of sinking air than a maximum 
convergence at 700 mb would. 

The existence of this trough be fore the development of the squall 
line would cause sinking and divergence at 700 mb and below and rising 
and convergence above 700 mb in front of the trough. Unfortunately, 
since the mesonetwork was not activated until 1700 GMT and the initial 
echo was recorded at 1530 GMT, whether or not such a state existed in 
the mesoscale environment is impossible to determine for this case. 

If the 700 mb trough is moved back in time at its average propogation 
from 1700 to 2000 GMT to where it would have been at 1530 GMT, how 
would it compare to the position of the initial echo? Fig. 21 shows 
this estimated position and the initial echo at 1530 GMT. It shows 
that the initial echo is almost exactly under the inflection point 
in front of the 700 mb trough. This is the area of upward motion 
above 700 mb with increasing vorticity and decreasing stability. 

This then strongly indicates that the storm was initiated aloft and 
worked down through the lower levels. 

Van Sickle (1969), in a study of the storms in the same area on 
28 May 1967, also came to the similar conclusion that the storms 
were initiated aloft. However, he found evidence of a mesoscale 
wave above 200 mb, therefore his initiation position was below and 
behind the wave. His wave was very similar in all aspects to the one 
found here, except that it was evident much higher in the atmosphere. 

2 . As a Possible Terminator of the Squall Line 

If the storm were initiated by upward motion above 700 mb, 
it would initially produce an echo with a relatively high base. As 
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FIG. 21, Estimated position of the 700 mb trough at 
1530 GMT and the initial echo, 
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the cloud grew it would definitely begin to interact with and reshape 
the mesoscale environment. As indicated by Takedo (1965), the 
building echo aloft would cause mass convergence upwind of it because 
of wind shear in the vertical. This then would cause an increase in 
the downward motion behind the cloud, and cause the base of the downward 
motion to rise above 700 mb. This in turn would cause an increased 
divergence below the inflection point near the surface. It has been 
shown that the squall line is moving faster than the 700 mb trough. 

As the echoes build, they not only cause convergence behind them aloft, 
but also are in turn controlled in movement by some average wind aloft, 
suggested in this case to be that near 500 mb. This was also the level 
of control found by Van Sickle (1969). As the cloud works up and down 
in the atmosphere and slowly moves away from its 700 mb ’’mother" (the 
inflection axis in front of the trough), it becomes the dominant 
feature in the area, its own motions slowly tending to dampen those 
caused by the environmental 700 mb trough. It soon becomes com- 
pletely dominant, setting up the large convergences near the surface 
where the necessary moisture to sustain it is located. During this 
stage the maximum rainfall is recorded. However, as the squall moves 
further away from the 700 mb trough, it soon moves into the area of 
the inflection axis immediately preceding the mesoscale ridge (about 
70 n.m. in front of the trough). Here, although definitely weakened 
by the squall itself, a weak divergence field remains. As the squall 
moves into it, the area of divergence acts as a destructive mechanism. 
Since the strongest convergence is near the surface in front of and 
near 550 mb behind the squall line, the 700 mb level of the squall seems 
the most susceptible to a divergence mechanism in the mesoscale 
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environment. Even a slight amount of forced divergence at 700 mb would 
cause a radical change in the squall. It would provide an outlet for 
the moist air from below, and cause the vertical motion pattern to 
decrease and actually reverse aloft. This would in return increase the 
divergence at 700 mb. The squall would in turn begin to ”die*' from aloft. 
Shortly the mesoscale trough-ridge system would once again be the dominant 
system in the area. 

The above proposed life cycle of the squall line as determined by 
the 700 mb mesoscale wave system admittedly is debatable and diffi- 
cult to substantiate because of the lack of data at the initiation 
stage and the movement of the squall outside the computational region 
as it begins to die at 2000 GMT. Also, knowledge of the amount and 
timing of the interaction of a squall line and its mesoscale environ- 
ment and their effects on one another is almost completely lacking. 
However, some evidence that the 700 mb mesoscale wave is the controlling 
mechanism can be found in the data. 

One indication that the squall was initiated aloft is the report 
of hail at 1600 GMT 20 n.m. southeast of LTS (directly in line of 
motion of the initial echo). This would indicate large vertical 
motions at and above the freezing level. The freezing level reported 
by LTS at 1700 GMT is between 600 and 650 mb. Yet comparison of the 
composite rainfall between 1615 and 1745 GMT (Fig. 10) to that between 
1745 and 1915 GMT (Fig. 13) shows considerably more rainfall in the 
latter period. This would indicate that an upward vertical motion 
was present above 700 mb at 1600 GMT but not necessarily below it, 
where the vast majority of precipitable water exists. 

Referring once again to the general fields calculated at all 
three times^ do they support the life cycle scheme presented? At 



84 



1700 GMT, tlie upper levels (Fig. 10) indicate the maximum vertical 
motion is oriented parallel to but slightly behind the squall line, 
almost directly over the inflection point at 700 mb. At the lower 
levels (Fig. 9) some divergence and sinking is still indicated below 
the center (youngest part) of the squall line near the inflection axis, 
yet reasonably large (5 to 10 x 10 ^ sec ^) convergences are noted in 
front of the squall. So, this would appear to be a growth stage of 
tlie squall, with the squall line and mesoscale environment having 
approximately the same magnitude of effect on the resultant 
ca Icula tions . 

By 1830 GMT, it appears the squall line is the dominant factor 
controlling the environment in the area. Fig. 13 indicates the 
maximum rainfall for all three periods considered is during this time. 
The upper level shows (Fig. 13) strong upward motions over the squall, 
but now moderate to strong (maximum of 60 mb hr ^) sinking exists 
between the trailing 700 mb trough and inflection axis. This supports 
the higher level of strongest convergence mentioned previously which 
was proposed to be the result of the blocking action of the great 
heights of the squall cells. (Another factor in favor of the largest 
cells at this time is the previously mentioned fact that the squall 
moved faster from 1830 to 2000 QIT than previously.) The lower 
levels (Fig. 12) show intensified divergence and sinking near the 
trough, and convergence and upward motion in the squall area, all 
indicating the squall line is in full swing. 

By 2000 GMT the squall line is past the 700 mb mesoscale ridge 
(Fig. 16) and definitely is weakened. The upper level sinking behind 
it has weakened sharply (maximum of 20 mb hr ^) and is now centered 
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behind the 700 mb trough. Except for a weak strip (5 mb hr of 
sinking across the center of the inflection axis preceding the trough, 
upward motion appears over the inflection axis. The convergence has 
decreased in both magnitude and area. The lower levels (Fig. 15) 
still show divergence and sinking under the trough and inflection axis. 
Yet, as mentioned previously, numerous cells are now developing through- 
out the area, and the rainfall isopleths appear to extend backwards 
into the area (Fig. 16), 

The above would seem to support the case that the mesoscale trough 
system is once more dominant and setting up similar fields to those 
supposed existing before the initial squall formed at 1530 GMT, 

No data was studied after 2000 OIT, However, radar pictures (not 
presented) indicated some organization of cells in a north-south line 
coming out of the area under the inflection axis at 2000 GMT. This 
organization was very weak however, and random cells remained through- 
out the area. 
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VI. CONCLUSIONS AND COMMENTS 



Analysis of the data for three times covering a 3-hour period of 
a squall line in central Oklahoma has shown that a mesoscale wave 
exists at 700 mb with a wavelength of 90 to 100 n.m. in the area at 
all three times behind the squall line. The wave is moving at nearly 
twice the speed of the 700 mb wind, but continually loses ground to 
the squall line, which is apparently being moved at about the speed 
of the 500 mb wind. 

The mesoscale wave causes convergence in front of it, with sink- 
ing motion at 700 mb and below, and rising motion above. Its extrapo- 
lated position at 1530 GMT (the time of the initial echo) indicates it 
is causing upward motion above 700 mb with increasing vorticity and 
decreasing stability over the area where the initial cell developed. 
This is proposed to be the cause of the squall line. 

As the squall line moves to a position approximately one-half of 
the 700 mb wavelength in front of its birthplace with respect to the 
700 mb wave it begins to die. Since the opposite conditions (sinking 
aloft and rising below 700 mb) exist here than those where the storm 
was initiated, it is felt the 700 mb pattern causes the death of the 
squall line. 

In between there appears to be a delicate interaction of the 
squall line and the mesoscale wave, a complete description of which 
is beyond the scope of this paper, and the science at this date. In 
general, the mesoscale wave appears to be the controlling factor at 
the initial and decay stages, and the mature squall line dominates 
the atmospheric fluctuations in between. 
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The squall line itself, at least during the time it dominates the 
tropospheric motions in its vicinity (mature stage), has a structure 
in general agreement with the classical picture of low level conver- 

-4 -1 

gence and upper level divergence. The strong convergences (10 sec ) 
in the lowest levels, where sufficient moisture to maintain the squall 
exists, are found in the vicinity of and preceding the squall. Upward 
vertical motion maxima are found in and near the strongest cells near 
450 mb. Divergence occurs above the jet stream, generally from 200 to 
150 mb. The strong convergences noted from 650 to about 400 mb seem 
to indicate that considerable entrainment may be taking place, although 
this aspect was not investigated. These large convergences aloft, 
generally behind the squall, cause strong sinking and associated low 
relative humidities in the areas immediately proceeding the squall 
line. They may indeed support the "blocking action" of well organized 
convection as postulated by Newton (1966) and Alberty (1969). These 
large convergences aloft, also noted by Fankhauser (1969) and Van 
Sickle (1969) for the case of 28 May 1967, merit more intensive study. 

Because the sounding network is always initiated at 1700 GMT, 
much of the possible valuable information that could be learned from 
the data gathered by it is lost, as on many occasions (including the 
one presented here) considerable activity has already developed by 
this time. Since a representative picture of the environment before 
the initial activity cannot be obtained, verification of the initial 
cause (in this case the 700 mb mesoscale wave) cannot be made satis- 
factorily. It seems that such an expensive and otherwise highly 
efficient network is wasting a valuable opportunity by not beginning 
observations prior to the initiation of convective activity. 
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APPENDIX A 



THE PROBLEMS OF ASCENSION RATE, RELEASE TIME, TIME LAG, 

AND BALLOON DISPLACEMENT 

As mentioned in Section III, the drift of a balloon from its launch 
position, the time it takes a balloon to reach a pressure level as com- 
pared to its launch time, varying launch times, and varying ascension 
rates of balloons (due to varying amounts of helium) are generally not 
considered problems on the synoptic scale because of the relatively 
small magnitude of these variations when compared to synoptic-scale 
values. However, on the mesoscale, because of the relatively short 
distances between stations (in this case 40-50 n.m.) and relatively 
short times between soundings (in this case 1-1/2 hours), these same 
variations are, relative to the mesoscale values, much larger and 
could possibly cause a misleading presentation of the atmosphere if 
not corrected for properly. Fankhauser (1969) has noted all these 
possible problems on the 28 May 1967 case and attempted to correct 
them whenever necessary by various means. 

In the case examined here, the problems of differing ascent rates 
and differing times of release of the balloons would appear negligible, 
considering the extreme care taken to meet the standard requirements 
set forth to be followed by this network, also noted by Fankhauser 
(1969). In order to determine whether or not these two possible 
problems were significant for this case, the release times and 
ascension rates of the soundings (25 for this 3-hour period) were 
carefully examined. Table I shows that the maximum variation from 
the predetermined release time was 15 minutes at TIK at 1700 GMT. 
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Table II shows that the average ascension rate to 150 mb for all releases 
was about 49 minutes, with a maximum variation for all 25 soundings of 
+ 9 minutes. It is felt that these facts show that varying times of 
release and ascension rates are generally negligible for this case. 

Now consider the possible problems of balloon displacement and time 
lag. Matsumoto and Ninomiya (1963) have noted that failure to take 
into account balloon displacement on the mesoscale can lead to errors 
of as much as 70^ in the computation of divergence in the upper levels. 
However, their argument is not considered too applicable here as they 
considered a three station Gaussian method of computing divergence, 
considerably different from the standard method used here and explained 
in Section IV. Nor did they mention the problem of time lag at the 
pressure levels. 

In order to determine if balloon displacement and time lag were 
problems in this case which could not be ignored, the following pro- 
cedure was followed. To construct a vertically consistent sounding 
at 1830 GMT, nine levels at 1700 GMT and 1830 GMT above 600 mb (550 
to 150 mb) were replotted, taking into account balloon displacement 
(i.e., plotting winds exactly where they were measured with respect 
to the reporting station). In general, rather light winds and/or 
winds that shift considerably with height occur at 600 mb and below. 
Tables II and III show that the errors associated with balloon dis- 
placement and time lag below 600 mb (average ascent time of 15 minutes 
to 600 mb and horizontal displacement of less than 7 km) are within 
the analysis limitations for a network of this density. These charts 
(550 to 150 mb) were analyzed and data taken using the 13 x 11 grid 
previously mentioned in Section I. 
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STATION 


1700 GMT 


1830 GMT 


2000 GMT 


SPS 


1651 


1828 


1952 


RIN 


1650 


1822 


MSG 


LTS 


1656 


1821 


1951 


FSI 


1653 


1837 


2009 


PVY 


1654 


1830 


1956 


COR 


1700 


1830 


2000 


CHK 


1651 


1819 


1950 


TIK 


1715 


MSG 


2000 


WAT 


1700 


1830 


2000 


STANDARD 


TIK 


CHK 


CHK 


DEVIATION 


-15 


+11 


+10 



TABLE I. Release times and maximum deviations 
from predetermined launch times for the rawin- 
sonde mesonetwork in central Oklahoma on 30 May 
1967. 
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LEVEL (mb) 


900 


750 


600 


450 


300 


150 


Average time to 
(min) 


2.21 


8.19 


15.1 


23.47 


34.68 


48.92 


Maximum time to 
(min) 


4.15 


10.66 


21.11 


30.12 


42.29 


58.44 


Minimum time to 
(min) 


1.32 


5.81 


10.95 


14.98 


27.85 


40.21 



TABLE II. Ascent times to various pressure levels of rawinsondes 
launched between 1700 and 2000 GMT on 30 May 1967 by the mesonet- 
work in central Oklahoma. 



LEVEL (mb) 


900 


750 


600 


450 


300 


150 


Average Displacement 
(km) 


1.35 


3.77 


6.88 


15.67 


28.93 


50.76 


Maximum Displacement 
(km) 


2.8 


9.4 


16.1 


00 

CM 


42.5 


68.8 


Minimum Displacement 
(km) 


0.4 


0.6 


2.2 


8.1 


19.4 


39.7 



TABLE III. Balloon displacement from launch site at various 
pressure levels of rawinsondes launched between 1700 and 2000 
GMT on 30 May 1967 by the mesonetwork in central Oklahoma. 
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Having balloon displacement solved, the time lag problem must also 
be solved, at least to the first approximation. Since the balloon 
displacement lias been taken into account previously, it is more correct 
now to linearly interpolate for time correction (rather than vice-versa). 
Again, Tables II and III show the time correction to be considerably 
more linear than balloon displacement. So, using the reported times 
that the several balloons reached each specific pressure level, an 
average time for each pressure level (for each of the two general 
sounding times) was computed. At 1830 GMT, it was determined that 
all the soundings were representative of the environment and usable. 

At 1700 GMT, because of the tardiness of the launch at Tinker (TIK) 
and an early release at Wichita Falls (SPS), these two times were 
deleted when determining the average times because of their large 
deviation from the average of the rest of the soundings. The rest 
of the soundings at both times at each individual level were within 
a maximum o f 14^ of the average time of the level, and most cases were 
within + 9 minutes of the average, or 10^. It was decided that the 14^ 
variations were within tolerance, generally in an attempt to maintain 
as much of the limited data as possible. The 1700 GMT TIK data was 
used as an approximation, since it appeared to fit the environment. 
However, little if any use could be made of the SPS sounding, since 
it was near or in, especially in the upper levels, the leading edge 
of the squall line. 

The general computer program, developed with advice from Dr. J, D, 
Mahlman, adjusted the u and v components of the wind field at 1830 GMT 
for time lag before being used for computa t ions . This was done by 
using the average times computed and stored for each level at both 
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1700 GMT and 1830 GMT as indicated below: 



r -\ ( 1830- 1 - ) 

B(1830) = B(t^) + [BCt^) - B(t^)J 

where B = any measured or derived variable, 

t^ = computed average time of level for 1700 GMT soundings, 

and t^ = computed average time of level for 1830 GMT soundings. 
The above adjustment is made in subroutine COMSND in the computer 
program . 

So, the u and v components of the atmospheric wind field at 1830 
GMT on 30 May 1967 in the grid area should be closely represented. 

Comparisons of the output fields using these adjusted u and v 
components and the unadjusted ones at 1830 GMT were made to deter- 
mine whether or not the data could be used in the general manner 
without adjustments for balloon displacement and time lag. Figs. 

22 - 25 show vertical profiles of adjusted divergence and vertical 
motion (omega) both uncorrected and corrected at four points in the 
grid area. For clarity, it is once more emphasized that all vertical 
motions presented in any figure in Appendix A are those determined 
by the adjusted kinematic technique (final step) as explained in 
Section IV-C. For these points, selected because they are in grid 
areas where significant values of vertical motion (+) are occurring, 
the uncorrected fields of u and v give a very close approximation to 
the profile as determined after corrections for balloon displacement 
and time lag. It is noted that magnitudes of both divergence and 
vertical motion vary somewhat, but the variability is small and in 
the same direction as the corrected profile, therefore causing no 
drastic change in the overall profile shape, certainly no more than 
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FIG. 22. Uncorrected divergence (xxx) and vertical motion 

(...) versus divergence ( ) and vertical motion ( ) 

corrected for balloon displacement and time lag at 1830 GMT 
on 30 May 1967 for grid point (5, 5) located in the Oklahoma 
rawinsonde mesonetwork. (Note: For location of grid points 

see Fig. 1). 
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FIG. 23. Uncorrected divergence (xxx) and vertical motion 

(...) versus divergence ( ) and vertical motion ( ) 

corrected for balloon displacement and time lag at 1830 GMT 
on 30 May 1967 for grid point (9, 7) located in the Oklahoma 
rawinsonde mesonetowrk. (Note: For location of grid points 

see Fig . 1) . 
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FIG. 24. Uncorrected divergence (xxx) and vertical motion 

) versus divergence ( ) and vertical motion ( ) 

corrected for balloon displacement and time lag at 1830 GMT 
on 30 May 1967 for grid point (10, 6) located in the Oklahoma 
rawinsonde mesonetwork. (Note: For location of grid points 

see Fig . 1) . 
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FIG. 25. Uncorrected divergence (xxx) and vertical motion 

(...) versus, divergence ( ) and vertical motion ( ) 

corrected for balloon displacement and time lag at 1830 GMT 
on 30 May 1967 for grid point (9, 5) located in the Oklahoma 
rawinsonde mesonetwork. (Note: For location of grid points 

see Fig . 1) . 
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the smootliing process does. It appears that fiad Matsiimoto and Ninomiya 
(1963) also taken into account tlie time lag of the balloon from obser- 
vation time along witli balloon displacement, instead of only displace- 
ment, considerably less error would possibly liave been noted in the 
upper Levels. Of the four profiles noted here, two (Figs. 22 and 25) 
show a change in sign of the vertical motion with height. It is 
interesting to note that the corrections caused practically no change 
in the level at which the change in sign occurred. Also, noting the 
divergence patterns, the same applies to all four profiles where there 
is a change in the sign of the divergence. The above seems to indicate 
that corrections for displacement and time lag are not necessary for 
this case. 

Leaving the vertical for a moment, examine the horizontal fields 
before and after correction. Figs. 26 and 27 show at 350 mb uncorrected 
and corrected omega fields at 1830 GMT respectively. A general look 
appears to indicate very similar patterns, except possibly in the 
southeast corner. Looking closely at grid points indicates that all 
maximum positive and negative omegas, although again varying slightly 
in magnitude, occur at the same grid points, except in two areas. 

These are the maximum negative center in the east central area and 
the maximum positive area in the southeast. The negative center has 
the same orientation on both figures and approximately the same 
magnitude. However, the maximum value is displaced one grid point. 

Since the actual maximum value could occur anywhere between the two 
grid points, one grid point variation is a definite limit of the 
system. Thus, the only problem area is the relatively large positive 
center in the southeast in Fig. 26 versus a smaller positive center 
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FIG. 26. The 350 mb omega field (mb hr ) calculated 
from the wind field which has not been corrected for 
balloon displacement and time lag at 1830 GMT on 30 
May 1967. 
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FIG. 27. The 350 mb omega field (mb hr ) calculated 
from the wind field which has been corrected for balloon 
displacement and time lag at 1830 GMT on 30 May 1967. 
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further to the southeast in Fig. 27. The error here can be represented 
by comparing the values from both figures at grid point (11, 4) which 
is just north of the positive center in Fig. 26 (see Fig. 1 for exact 
position of any grid point). In Fig. 26 it is +35, versus -10 in Fig. 
27. Therefore, there is a definite need to correct for balloon dis- 
placement and time lag here. Why? 

Pose the question of where one would expect the time lag not to 
cancel out the balloon displacement. It would seem that since the 
time lag varies almost linearly with height it would fail to approxi- 
mately balance the balloon displacement when the balloon displacement 
increased rapidly and non- linearly over the same height. Where would 
an excessive balloon displacement occur? It would occur when there 
was a large vertical shear or a considerable change in wind direction 
with height. The former would cause a non-linear increase in displace- 
ment, and the latter a change in displacement direction which would 
make the fields measured downstream not compatible with that which 
occurred at the earlier time at the position of interest. Do such 
conditions occur at time 1830 GMT? If they do, is there any connection 
with the problematic southeast area noted earlier? 

An examination of the isotach and streamline fields at 1830 GMT 
above 600 mb shows several positions where considerable change in 
speed or direction with height occurs. Since errors are additive 
with height and the scheme used for determining omega (Section IV-C) 
causes more correction with height because of radiosonde errors, 
concentrate on the first few levels above 600 mb where changes occur. 

Figs. 28 and 29 show the streamline and isotach fields at 1830 GMT 
for 500 and 400 mb respectively. Comparison of the two shows several 
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FIG. 28. Streamlines ( — ^>) and isotachs ( ) (m sec ) 

for 500 mb at 1830 GMT on 30 May 1967. 
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FIG. 29. Streamlines ( — ^>) and isotachs ( ) (m sec 

for 400 mb at 1830 GMT on 30 May 1967. 
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changes between the two levels. In particular, the wind backs in the 
northwest, decreases in the east central area, and veers and decreases 
in the southeastern area with height. It was noted previously that the 
southeast area was the problem area, as noted by Figs. 26 and 27. The 
northwest and east central area were considered acceptable. Glancing 
again at Figs. 26 and 27, pick three points ((2, 8), (11, 5), (11, 4)) 
which correspond with the northwest, east central, and southeast areas 
noted above. Figs. 26 and 27 show a general agreement in omega for 
(2, 8) and (11, 5), but strong disagreement as noted previously at 
(11, 4). Table IV gives wind directions and speeds for these three 
grid points at 500 and 400 mb. Table IV and Figs. 26 and 27 then 
indicate that a change of 45° and 7 m sec ^ at (II, 4) per 100 mb is 
significant, but changes of 15° and 12 m sec ^ at (11, 5) and 30° and 
1 m sec ^ at (2, 8) per 100 mb are not. Vertical plots of these three 
grid points above 600 mb are shown in Figs. 30 (2, 8), 31 (11, 5), and 
32 (11, 4), Comparisons of vertical motions and divergence before and 
after corrections for displacement and time lag in these figures 
indicate a generally good agreement in Figs. 30 and 31, but Fig. 32 
indicates that corrections must be made to get a representative 
picture of the environment. Therefore, Table IV, Figs, 26 and 27, 
and Figs. 30, 31, and 32 are all in agreement. 

The above indicates that changes in speed of less than 15 m sec ^ 
or direction of 30° or less per 100 mb by themselves do not necessi- 
tate making corrections for balloon displacement and time lag. Seldom, 
at altitudes above 600 mb will changes greater than these occur, except 
possibly near the tropopause where necessary corrections are so great 
because of radiosonde error that results are not particularly dependable 
a ny wa y . 
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FIG, 30. Uncorrected divergence (xxx) and vertical motion (...) 

versus divergence ( ) and vertical motion ( ) corrected for 

balloon displacement and time lag above 600 mb at 1830 GMT on 
30 May 1967 for grid point (2, 8) located in the Oklahoma 
rawinsonde mesonetwork. (Note: For location of grid points 

see Fig. 1). 



107 




FIG. 31. Uncorrected divergence (xxx) and vertical motion 

(...) versus divergence ( ) and vertical motion ( ) 

corrected for balloon displacement and time lag above 600 
mb at 1830 GMT on 30 May 1967 for grid point (11, 5) located 
in the Oklahoma rawinsonde mesonetwork. (Note: For location 

of grid points see Fig. 1). 
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FIG. 32. Uncorrected divergence (xxx) and vertical motion 

(...) versus divergence ( ) and vertical motion ( ) 

corrected for balloon displacement and time lag above 600 mb 
at 1830 GMT on 30 May 1967 for grid point (11, 4) located in 
the Oklahoma rawinsonde mesonetwork. (Note: For location of 

grid points see Fig. 1). 
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A final check of the 1830 GMT data was made to get a comparison of 
the u and v components of the uncorrected data versus that corrected 
for balloon displacement and time lag. This check was a Root Mean 
Square (RMS) error analysis for several levels carried out by computer. 
Table V shows the results of this check, showing the RMSu and RMSv 
errors for several levels and the average u and v components for the 
uncorrected analyses at a few levels. The RMS values are all re la- 
tively small. The worst case shown is the RMSu at 300 mb as compared 
to the average u value, and that is about 30^, but the RMSv at 300 mb 
as compared to the average v value is of the order of 15^. As a rule, 
the RMS values are less than 20^ of their average component values at 
all levels. Considering the density of data (9 stations) and the 
possible 10 to 20^ source of error in analysis techniques alone, the 
results from Table V would indicate that no gross errors are made by 
not correcting for balloon displacement and time lag at 1830 GMT. 

Therefore, all of the above checks indicate that only a combination 
of direction and speed changes of greater than approximately 35° and 
10 m sec ^ per 100 mb, or individual changes of greater than 15 to 
20 m sec ^ in speed or 40 to 45° in direction per 100 mb necessitate 
a correction for balloon displacement and time lag. These are rather 
severe changes, and may be unusual, as is the case of the mesoscale 
data used here. 

The possible erroneous calculations have been noted at 1830 GMT, 
and the corrected values have been used in this study. The streamline 
and isotach fields at 1700 and 2000 GMT were checked for possible 
problem areas using the above stated criteria. Except for a few 
possible problem areas near 200 mb (the jet stream maximum problem 
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at specific pressure 



already noted), both 1700 and 2000 GMT data appear to meet the criteria 
limitations above. Therefore, no correction for balloon displacement 
and time lag was deemed necessary. 



If the criteria stated here can be considered a first approximation 
for a general rule, then hopefully, since this case appears to be one 
of reasonably large wind speeds aloft, a more general case might fit 
the criteria even better, and show even less error in computations 
when no corrections for balloon displacement and time lag are made. 

The advantages of being able to use data without corrections are easily 
seen. Some of these are : 

(1) considerable savings in time; 

(2) ability to use all data periods reported (correction scheme 
causes loss of at least one time period due to interpolation require- 
ments ) ; 



(3) and one is not dependent on the assumption of linear changes 
between soundings to construct a new sounding. 

Because of our limited ability to measure mesoscale phenomena or 
to separate mesoscale and microscale motions and their possible inter- 
actions, it seems hardly worthwhile and possibly misleading to use 
detailed techniques to correct admittedly crude data, if the data 
meets certain general criteria and therefore can be used as is. 
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APPENDIX B 



RESULTS OF THE COMPLETE VORTICITY METHOD (VOROMG) 

The complete vorticity method, labelled VOROMG in the computer pro- 
gram, was used in an attempt to determine vertical motion fields at 

dc 

r 

1830 GMT (one time was all that was possible because of the term - 

found in FUNl) . 

The same program, with minor alterations, was used by Dr . J. D, 
Mahlman to determine omega fields on the synoptic scale. He obtained 
total convergence in less than five cycles. The fields of omega 
obtained from this method (VOROMG) compared favorably with fields 
obtained by a separate omega equation method and those determined 
by the adjusted kinematic method (KINOMG) . However, the VOROMG omegas 
were somewhat larger on the average than those determined by KINOMG. 
Mahlman also found that the twisting term (FUN2) , which is shown to 
be negligible on the synoptic scale by scale analysis, is in reality 
quite often a significant contribution to the final answer. This is 
in good agreement with the previously noted results of Yanai and Nitta 
(1967). 

The success of the vorticity method on the synoptic scale was 
encouraging. However, all attempts to obtain convergence on the 
mesoscale failed. The method always diverged on the mesoscale, even 
though a first approximation using only FUNl (equation without the 
twisting term) showed that the requirement that omega be zero at the 
top of the atmosphere (150 mb) had been met. This same first approxi- 
mation indicated enormous (500 to 1000 mb hr ^) vertical motions in 
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the middle troposphere, with a few cases of adjacent values being near 
this size and of opposite signi The problem apparently was caused by 
the advection of vorticity in FUNl. Printouts of vorticity, especially 
in the lower and middle troposphere, showed much more complex fields 
than was the case for divergence. Apparently, cases of adjacent large 
values of vorticity of opposite sign, when acted upon by the vorticity 
advection term, caused fatally large values. These extremely large 
and sometimes random values of omega shown by the first approximation 
frustrated any possible attempt towards convergence. Attempts to smooth 
the vorticity field and to control the size of absolute vorticity (Q 

a 

found in the denominator of both FUNl and FUN2) were not successful. 

Even though the method failed on the mesoscale, its use on the 
synoptic scale may prove useful and very informative, especially in 
regard to the size of the formerly disregarded twisting term. 
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OPTERMIMjSTICN gf vorticity, oivergfncf, and vertical mGTIGM by thf 
ADJUSTED KINEMATIC AND COMPLETE VORTICITY mcthopS. 
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If CONTINUE 



TAKE OUT IF U8AR AND VBAR NOT WANTED 
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NVERT TO smaller FIELD AND PRINT VORTICITIES. THIS MUST BE 
NE TO ALL FIELDS TO MEET THE REQUIREMENTS OF SUBROUTINE 
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